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By studying periodic points for rational maps on d with p invariants, we show that they form 
an invariant variety of dimension p if the periodicity conditions are 'fully correlated', and a set 
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, of isolated points if the conditions are 'uncorrelated'. We present many examples of the invariant 
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varieties in the case of integrable maps. Moreover we prove that an invariant variety and a set of 



' isolated points do not exist in one map simultaneously. 
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When a map is given there is no immediate way to foresee the fate of its iterations or whether it is integrable 
or nonintcgrable. If the map is integrable the effect of a small difference of initial values remains small, while it 
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O |. behaves chaotically for some initial values if the map is nonintegrable. In spite of this large difference one cannot 
' easily distinguish the two cases. 

a- 

It would be desirable to know a way to distinguish integrable maps from nonintegrable ones just by investigation 
of the first few steps. The purpose of this article is to present a step of finding such a criterion in the case of higher 
, dimensional maps. 

5-H ' .„. 

In the theory of ordinary 2nd order differential equations, nonexistence of branch point singularities, which 
depend on initial values, indicates integrability of the equation (the Painleve test). The method of singularity 
confinement was proposedQ to replace this test and apply to discrete maps, but it does not always work[i]]. As far 
as one dimensional maps are concerned there have been intensive studies mainly from the mathematical point of 
view^, Q |{]]. It is known, for example, that the appearance of a Julia set, the closure of the set of repulsive periodic 
points, characterizes nonintegrability of the map. On the other hand very little is known about higher dimensional 
maps. In practice the best we can do at present is to study statistical information of the mapsj" 
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We are interested in the transition of a map from a nonintegrable one to integrable one when a parameter in 
the map is changed continuously. In this paper we pay attention to the behaviour of periodic points of higher 
dimensional rational maps. If the map is nonintegrable we shall find a set of isolated points with fractal structure 
as a higher dimensional counterpart of the Julia set. We would like to know what object appears when the map 
becomes integrable by adjusting the parameter. 

Before we present the results of this paper let us prepare some words. We assume that the map has p (> 0) 
invariants. The periodicity conditions are called 'uncorrelated' if they determine all positions of the periodic points 
dependent on the values of the invariants. They are called 'fully correlated' if they do not determine the position of 
the periodic points but impose relations on the values of the invariants. An 'invariant variety of periodic points' is a 
variety determined by the invariants of the map alone, such that every point on the variety can be an initial point of 
the iteration, which stays on the variety before it returns to the same point. We would like to show in this paper the 
following results: 

• Proposition 1 

A set of periodic points form a set of isolated points if the periodicity conditions are uncorrelated, and an 
invariant variety of dimension p if they are fully correlated. 

• Proposition 2 

An invariant variety of periodic points and a set of isolated periodic points do not exist in one map simultaneously. 

In order to clarify our argument we study in detail the three dimensional Lotka-Voltcrra (3dLV) map in §2 as 
an example. We will see that the periodicity conditions, which arc 'uncorrelated' in generic systems, become 'fully 
correlated' in the integrable maps and one gets full varieties of periodic points, instead of isolated points. The proof 
of our propositions will be given in §3. We present in §4 more examples of invariant varieties of periodic points of 
higher dimensional Lotka-Volterra maps and Toda lattice maps. When the dimension d is larger than 5 the invariant 
varieties become intersections of d — p algebraic varieties. Some examples will be given explicitly in the cases of d = 5 
and 6. 

The existence of an invariant variety of some period does not guarantee integrability of the map, but our proposition 
2 tells us that there is no set of isolated periodic points. On the other hand, if there exists a set of isolated periodic 
points of some period, the map can not have an invariant variety but will have infinitely many sets of isolated 
periodic points, a typical phenomenon of nonintegrable maps. Further discussion about the criterion, however, will 
be impossible unless we know a precise notion of integrability of a map. 

Although our results in §3 hold in general rational maps, we presented mainly those examples of integrable maps in 
§2 and 4 which belong to the series of Lotka-Volterra maps. The integrability of this series is guaranteed in the sense 
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that all solutions are given explicitly in terms of the r function of the KP theory. They are not all of the integrable 
maps but include a large part of known integrable maps. We explain in the Appendix the relations of the series of 
Lotka-Volterra maps, the Toda lattice and the r function of the KP theory. There is no reason not to study other 
types of integrable maps. 



II. A STUDY OF 3 DIMENSIONAL LOTKA-VOLTERRA MAP 

Consider a rational map on C d 

x=(x u x 2) ...,x d ) -> X = (X u X 2 ,...,X d ) =:I (1) , (1) 

We would like to investigate the behaviour of the sequence: x — > XW -> J( 2 ) -> .... In particular we are interested 
in the transitions between integrable maps and nonintegrable ones. Many discrete maps have been known to be 
completely integrable|8i 13. lift ll lL Il2j | . Starting from such an integrable map we deform it by introducing a continuous 
parameter. The map, after the deformation, will be nonintegrable in general. We compare these two systems by 
changing the parameter and studying their differences. 
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For an illustration let us study the following example|l3( with (xi, x 2 , X3) = (x,y,z) and a parameter a G C. 

1-y + yz 1 - z + zx 1 - x + xy 
X = x- ■ ha, Y = y- ■ h 2a, Z = z- ■ . 2 

1 — z + zx 1 — x + xy 1 — y + yz 

The number of independent terms including x in the numerator of X is four. If we repeat the map it becomes 41 in 
the numerator of X^ 2 \ If we continue repeating it we get 734 of such terms in the numerators of X^ and more that 
20,000 ofXC' etc.. 



In order to measure the complexity of the map we can use the algebraic entropy[2j,|g. Here we define it by 

E := lim — \nD n 

n — >oo Ti 

with D n being the highest degree of powers of the initial variable x in the numerator after the nth map. In the case 
of J5J, E = In 3. It is known that the map behaves chaotically when E is non-zero in the sense that we fail to predict 
the long-term behaviour for some initial values. 

The situation changes significantly when a = in J5J. In this particular case the map (J2J) is called the discrete 
Lotka-Volterra map, and is known to be completely integrable|ll|. In fact all solutions are given explicitly in terms 
of a r function of the KP hierarchy for arbitrary initial values, as shown in the Appendix. 

Let us see what happens in the map J5J) when a = 0. If we repeat the map the number of independent terms 
including x in the numerator of X^ is 10. If we continue repeating it we get 68 and 300 such terms in the numerators 
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of X^ and X^\ instead of 734 and the number of order of 2 x 10 4 of the case of a =^ 0, respectively. This large 
difference owes to the cancellation of factors in the numerator and the denominator in each step of the map when 
a = 0. Correspondingly the highest power of x in the numerator of X^ n ' changes slowly as 1, 3, 7, 11, 17, • • •, and 
hence the entropy is zero. 

We introduced the parameter a as simply as possible to emphasize that a small deformation from the integrable 
regime may cause a large difference. Needless to say many other deformations could be chosen. Only if we are very 
careful do we have a chance to generalize the integrable Lotka-Volterra map without introducing chaos. 



The algebraic entropy is correlated with an increase of the number of periodic points of the map. In a one dimensional 
map the Julia set, which characterizes nonintegrability of the map, is defined as the closure of the set of repulsive 
periodic points. Therefore it will be quite natural to study the behaviour of periodic points to see an indication of 
the integrability of the map. 

If the map returns to the initial point after n steps visiting other points only once, we call the point a point of 
period n. A point of period one is called a fixed point. A periodic point is called attractive (repulsive) if a small 
neighborhood is mapped closer to (further from) the point. Otherwise the periodic point is called neutraljj]]. 

To see the problem closer let us first calculate the fixed points of the map @ when a ^ 0. There are two fixed 
points on the (x,y) plane in C 3 : 



. , / 2a + 3 + V4a 2 + 4a + 9 2a + 3 - \/4a 2 + 4a + 9 \ 

[x,y,z) = , , , (3) 



/ 2a + 3 - V4a 2 + 4a + 9 2a + 3 + sy/4a 2 + 4a + 9 \ 
{x,y,z) = I , , 1. (4) 

If we repeat the map we obtain 288 isolated points of period 2. In general the number of solutions to the set of 
periodicity conditions 

iW = x, yC") = y, ZW = z (5) 

increases very rapidly as n increases. After every step of the map we will get a new set of finite number of 'isolated 
points' in the space of initial values x — (x, y, z). From the construction of this set of points, it has a fractal structure 
as we repeat the map. If we continue this procedure infinitely many times the closure of the set of points might 
become dense. 

Now let us turn to the case of an integrable map. In the case of a = in © we find 'lines' of fixed points. In fact 
all points on the four lines 

{x = y = 0} U {y = z = 0} U {z = x = 0} U {x = y = z] (6) 
satisfy the conditions (X {1 \Y {1 \Z^) = (x,y,z). This surprises us because we expect only two points at 
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(x, y, z) = (3/2, 0, 0), (0, 3/2, 0) in the limit of a -> of © and Q . 



If we repeat the map when a = 0, we find, instead of 288 points, a hyperbolic surface in C 3 defined by 

{.s + l = 0}. 



(7) 



Here 



s:=(l-x)(l-y)(l-z). 



(8) 



is one of two invariants of the 3dLV map 

This phenomenon continues to the cases of higher period. Namely, if we denote by r the other invariant of the map 



r := xyz, 



(9) 



we find 



{r 2 + s 2 - 



rs + r + s - 



1 = 0} 



(10) 



for the period 3 case, 



{r 3 s + s 3 ~ 3rs 2 + 6r 2 s + 3rs - r 3 + s = 0} 



(11) 



for the period 4 case and 



{rV - r 3 s 2 - 6r 4 s 5 + 10r 3 s 6 + 3s 5 r + s e + s 5 + 3rV - 3r 5 s 3 - 6r 4 s 3 
-r 6 s 3 + 3r 5 s 4 + s 4 + 21s 4 r 2 + 6s 4 r + r 3 s 7 + s 7 + 27s 5 r 2 - 3s 6 r - r 3 s 5 
+21r 2 s 6 - 10r 3 s 3 - 6rs 7 + s 8 = 0}, 



(12) 



for the period 5 case, respectively. It is remarkable that the surfaces are characterized only by invariants of the map. 
These are the invariant varieties of periodic points, since any point on the surfaces can be an initial point and the 
following iteration of the map remain on the same surface before it comes back to the initial point. It will be possible 
to find more surfaces associated with higher periods if it is necessary. 



Now we are going to study the nature of the invariant varieties of periodic points just we have found. The behaviour 
of a neighbourhood of a periodic point is determined by multipliers Ai, A2, A^, i.e., the eigenvalues of the Jacobi 
matrix of the mapl^]]. 



/ dXi gXi _ _ _ dXi v 

dx\ dx2 dx n ' 



J 



0X 2 dX 2 
dxi dx2 



dX 2 

dx n 



(13) 



In the directions with |Aj| < 1 the map is stable while it is unstable in the directions with |Aj| > 1. If |Aj| = 1, the 
point is said to be neutral along the jth direction. 

If we calculate eigenvalues for the map J5J with o^Owe will see no particular regularity since they are governed 
by the value of a. To study the case of a = we notice that three equations of the periodicity condition J^J must be 
satisfied by a single constraint on the (x, y, z) for the periodic points to form an invariant surface. In other words an 
invariant variety of periodic points of period n emerges when the nth image of the map is given as 

Xjf>(x) = Xj + Uj {x) ln {x) j = 1, 2, 3, (14) 

with Uj's being regular functions when j n = 0. After some manipulation we can convince ourselves that the examples 
0, ijTUI) . ijTT|) and (|T2)) of the 3dLV correspond to 72, 73, 74 and 75, respectively. 
When l|14|l holds, the multipliers at a point on the variety v((j n )) are given by solving 

(A-l) 2 (A-det J (n) ) =0, 

with 

3=1 3 

Therefore it is guaranteed that every point on the variety is neutral along the variety as it should be expected. 
Moreover J^ n ' = 1 holds in the case of the 3dLV. Therefore all points on the invariant surfaces are neutral in all 
directions. 

In addition to (|10J) there are lines, in the case of period 3, specified by 

{x = y=l} U {y = z = l} U {z = x = 1}, (15) 

and 

{x = y = —oj} U {y = z = —u} U {z = x = —lu} (16) 

U {x = y = -Lj 2 } U {y = z = -uj 2 } U {z = x = -w 2 }, (17) 

where lu is the cubic root of 1. These lines are, however, peculiar in period 3 and not to be considered in the same 
level of the surface specified by (|10fl . In H15fl the union of the three lines should be conserved. Namely the map does 
not stay on one line, but goes around following the rule 

(1, 1, t) - (t, 1, 1) -> (1, t, 1) - (1, 1, t) -» • • • (18) 

with t € C. This is due to the discrete symmetry special to the 3dLV mapjl^. A point on one of the six lines given 
by i|16fl and (|17fl will not be mapped to the same line either. It corresponds to one of the three singularities at 

x(l-y) = l, y(l-z) = l, z(l-x) = l 
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of the map. If we map one of the points, say (t, 1 — 1/t, z), on the surface x(l — y) = 1 we have a sequence 

t,l - j,z\ -> (l,oo,0) -> (-oo, 1,0) ^ (l-j,t,z 

Repeating the map 6 times the point returns to the initial point. It has period 3 if we restrict to y = x, i.e., when 
t = —uj or — uj 2 , which are nothing but cases of (fTfi)) and fTTjl . 



III. THE NATURE OF PERIODICITY CONDITIONS 

We study in this section the nature of periodicity conditions of a rational map on C d which has p (> 0) invariants. 
The periodic points of period n will be found by solving 

X^= Xj , j = l,2,...,d. (19) 

If Hi(x), Hi{x), ■ ■ ■ , H p (x) are the invariants, the solutions of (|19|) are constrained on an algebraic variety of dimension 
d — p specified by the set of equations 

H i (x) = h i , i=l,2,...,p. (20) 

Here h\,hi, ...,h p are the values of the invariants determined by an initial point of the map. Let us denote this variety 
by V(h), i.e., 



V(h) = [x\ H i (x) = h i , i = l,2,..., p y 



The problem of finding periodic points is equivalent to finding an ideal generated by the set of d + p functions 
{Xj (x) — Xj, Hi(x) — hi}. Since the existence of the invariants enables us to eliminate p components of x from 
(|19fl . the ideal reduces to the pth elimination ideal generated by certain functions satisfying 

r ( °' ) (h 1 ,h 2 ,...,h p ,y 1 ,y2,...,yd-p) = 0, a = 1,2, ...,d-p. (21) 

Here by y\,yi, ...,yd-p we denote the variables which parameterize the variety V(h) after the elimination of the 
components of x. 

For an arbitrary set of values of hi, h p , the functions T^*\h, y) define an affine variety, which we denote 

V((T n )), i.e., 



V((T n )) = {y| rW(h,y) = 0, a= 1,2, ...,d - p}. 



In general this variety consists of a finite number of points, hence zero dimension, corresponding to the solutions 
to the d — p algebraic equations pifl for the d — p variables yi,y2, ■■■,yd-p- Once the values of 1/1,2/2, ■■■,yd~p are 
decided by solving l|21|l , the location of a periodic point on C d will be determined from the information of the values 
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of hi, hi, h p . In this way we obtain a number of isolated periodic points of period n. This is the precise meaning 
of the 'uncorrelated' periodicity conditions which we introduced in §1. Needless to say this case includes a map with 
no invariant. 

There are possibilities that the equations = impose relations on hi, hi, • ••> h p instead of fixing all y a 's. If 
s is the number of independent constraints imposed on h^s, the same number of y a 's are left free. Accordingly the 
dimension of the affinc variety V((L n )) becomes s. Every sequence of the map starting from a point on V((T n )) will 
return to the point after n steps. 

When s = d — p, all y a 's are free, while the values of the invariants are constrained. Let us denote by V n (h) 
the variety V(h) whose arguments (hi, hi, ■ h p ) are constrained by the d — p relations. Then the affine variety 
V((T n }) coincides with the variety V n {h) itself. In other words every point on V n (h) is a periodic point of period 
n. Thus all other maps are forbidden on V n (h), in a strong contrast to generic cases where s < d — p and none of 
map of other type is excluded on it. In this particular case, we denote by jn (hi, hi, h p ) the d — p functions 
r^* (hi, hi, hp, yi, 2/2> •••) Vd-v) to emphasize independence from the variables yi, yi, ya-p- The functions 
impose the conditions on V(h), but do not specify a point on it. 

Let us further replace hi by Hi(x) in 7„ and write the periodicity conditions as 



The expression (1221 enables us to consider the constraints on the invariants as constraints on the variables x. Un- 
der this new interpretation of the constraints the values of hi, hi, h p are not fixed but only relations among 
Hi(x), Hi(x), Hp(x) are imposed to decide the periodic points. If h^s fulfil the relations jn (hi, hi, h p ) = 0, 
all points x satisfying 1221) are periodic points of period n. We denote by v((7n)) the affine variety generated by the 
functions 7^ (Hi(x), Hi(x), H p (x)), and distinguish it from V((T n )). Namely we define 



Note that w((7n)) is a subvariety of C d , whereas V((T n )) is a subvariety of V(h), since the values of hi, hi, h p are 
fixed in the latter. 

The significance of defining v((-f n )) lies on the fact that for a point to belong to v((j n )) is a sufficient condition for 
the point being a periodic point of period n. Every point on f((7 n )) can be an initial point of the periodic map of 
period n. This is true only in the case of s = d — p. If s is less than d — p, the constraints on the invariants are not 
sufficient to decide periodic points. We can summarize properties of w((7n)) as follows. 

1. The dimension of v((^ n )) is p. 

2. Every point on v({^ n )) can be an initial point of the periodic map of period n. 



7, 



i a \H 1 (x),Hi(x),...,Hp(x)) = 



a = 1,2 



, ...,d-p. 



(22) 
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3. All images of the periodic map started from a point of u((7n)) remain on it. 

4. w((7„)) is determined by the invariants of the map alone. 

We have already shown many examples of v((j n )) in the case of the 3dLV map. We called v((j n }) in §1 an 'invariant 
variety of periodic points'. Since the 3dLV has 2 invariants, the number of independent constraints s = d — p is one, 
hence the dimension of the invariant varieties is 2. We shall present more examples in §4 including the cases of 
d-p = 2. 

By studying the periodicity conditions 119(1 we have shown that the nature of periodic map changes largely 
dependent on the constraints imposed on the invariants. In generic case, where values of the invariants can be chosen 
arbitrary, the periodic points appear isolated on C , We have called in §1 the periodicity conditions being 'uncorre- 
cted' in this case. When the invariants are constrained by d — p conditions, on the other hand, the periodic points 
form an invariant variety of periodic points. This type of conditions was called 'fully correlated'. We would say the pe- 
riodicity conditions being 'correlated' if they are not uncorrelated. The results were summarized in Proposition 1 of §1. 

We now suppose that, for a map with p invariants, the periodicity conditions of period k is fully correlated, 
hence there exists an invariant variety f((7fc)) of periodic points. At the same time we assume that the periodicity 
conditions of period n (^ k) are uncorrelated, i.e., dependent on the variables yi,y2, ••• J 2/d— p- All of the solutions 
to (hi, /12, ...,h p ,yi,y2, Vd-p) — determine points of period n for a set of values of h±, hi, h p . It must be 
true even when h\,h-2, ...,h v are chosen to satisfy ^^\hi, h%, h p ) — for all a. Recall, however, that all points on 
u((7fe)) are points of period k irrespective of y\,yi, ■■■,Vd-p- This contradicts to our assumptions that n ^ k and the 
periodicity conditions of period n are uncorrelated. Therefore the following statements are true for a rational map on 
C d . 

• If there exists a set of fully correlated periodicity conditions, say of period k, all other periodicity conditions of 
period n =/= k are correlated. 

• If there exists a set of uncorrelated periodicity conditions, there exists no set of fully correlated periodicity 
conditions. 

The second proposition of §1 follows after these results immediately. 

IV. INVARIANT VARIETIES OF HIGHER DIMENSIONAL MAPS 

Based on the study of the 3 dimensional Lotka-Volterra map in §2, we obtained in §3 some results which can be 
applied to all rational maps in higher dimensions. Our propositions do not guarantee that an existence of an invariant 
variety of some period implies integrability of the map. We have not proven either that every integrable map has an 
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invariant variety of periodic points. Nevertheless it will be worthwhile to study many known integrable maps and see 
if they have invariant varieties of periodic points. 

We present in this section some examples of invariant varieties derived from integrable maps whose dimensions are 
higher than 3. In all examples we shall see that the dimension of the invariant variety coincides with the number p 
of invariants of the map, in agreement with our result of section 3. The first examples are the Lotka-Volterra map 
and the Painleve V map whose dimension is 4. The dimension of the invariant varieties is 3 in both cases, while 
the number of the invariants is also 3. The second example is the 5 dimensional Lotka-Volterra map which has 3 
invariants. We shall find an invariant variety of dimension 3 as we expect. The final example is the 3 point Toda 
lattice map, which is equivalent to the 6 dimensional Lotka-Volterra map. The dimension of the invariant variety is 4, 
in agreement with the number of the invariants. As the dimension of the map and/or the degree of period increases, 
the procedure of finding the invariant varieties becomes harder, since we must manipulate algebraic formulae of many 
variables and/or of higher degree. 



A. 4 dim Lotka-Volterra map 



The 4 dimensional Lotka-Volterra map is defined by 



Xi = xi 



X 3 = x 3 



1 - X 2 - X 3 + X 2 X 3 + X 3 X4 
1 — X 3 — X4 + X3X4 + X4X1 
1 — X4 — X\ + X4X1 + X1X2 



X 2 = x 2 



X4 = X4 



1 — X 3 — X4 + X 3 X4 + X4X1 



1 — X4 — x\ + X4X1 + x\x 2 ' 

1 - X\ — X 2 + X\X 2 + X 2 X 3 
1 - X 2 - X 3 + X 2 X 3 + X 3 X4 ' 



1 - Xi — X 2 + XiX 2 + X 2 X 3 



There are three invariants of the map given by 



r = 



X\X 2 X 3 X4, 



t=(l 



xi - x 3 )(l - x 2 - Xi) 



U = X 2 X 3 X4 + X 3 X4X\ + X4X\X 2 + X\X 2 X 3 — X\X 3 — X 2 X4- 



We find an invariant variety of periodic points 



v(( l2 )) = {x\ t + 1 = 0} 



for the period 2 case, and 



w«73» = {x\ t 3 +t 2 r -t 2 u + t 2 -u 2 + 2rt - u + r + t = 0} 



for the period 3 case. 
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B. Painleve V map 



As an example which does not belong to the Lotka-Volterra series let us examine a discrete analogue of the Painleve 
V equation given by^J 



Xi := xi 
X 3 := x 3 



1 - X2 + X 2 X 3 - X2X3X4 
1 — X4 + X4X1 — X4X1X2 
1 — X4 + X4X1 — X4X1X2 



' 1 - X2 + x 2 x 3 - X 2 X 3 X 4 
Using the invariants of this map 



X 2 := x 2 
X4 := X4 



1 — X 3 + X3X4 — X3X4X1 
1 — X\ + X\X 2 — XiX 2 X 3 

1 — x\ + £1X2 — xix 2 x 3 
1 — X3 + 0:32:4 — 2:32:42:1 



r = 2:1X2X3X4, s = (1 - xi)(l - X2)(l - X3)(l - X4), v = (1 - x 2 x 4 )(l - X1X3), 
we find an invariant variety of periodic points 

"((72)) = {x\ s + v = 0} 

for the period 2 case and 

v((l3)) = {x\ (s + v) 2 -s(l-r) 2 = 0}. 

for the period 3 case. 



(23) 



(24) 



C. 5 dim Lotka-Volterra map 

The Lotka-Volterra map x — > X of dimension e? is defined by solving 

X j (l-X j _ 1 )=x j (l-x j+1 ), j = l,2,...,d (25) 

for X — {Xi,X 2l ...,Xd) under the condition Xj + d — Xj. We show in the Appendix that the number of invariants of 
this map is 

V = [d/2] + 1, 

where [d/2] is d/2 if d is even and [d/2] = (d — l)/2 if d is odd. The cases of d = 3,4, which we studied already, cor- 
respond to p = d—1. The problem of finding invariant varieties of periodic points becomes more difficult as d increases. 



By solving 125|) in the case of d = 5 we obtain 



1 


- x 2 


- x 3 


— X4 - 


V x 2 x 3 - 


h X2X4 - 


h X3X4 - 


h x 4 x 5 


- X2X3X4 


- X2X4X5 - 


1- X 2 X3X 4 X 5 


1 


- x 3 


— X4 


- x 5 - 


1- X3X4 + X3X5 - 


(- X4X5 - 


1- X5X1 


- X3X4X5 


- X3X5X1 - 


1- X3X4X5X1 


1 


- x 3 


— X4 


- x 5 - 


1- X3X4 - 


h X3X5 - 


F- X4X5 - 


1- X5X1 


- X3X4X5 


- X3X5X1 - 


f X3X4X5X1 


1 


— X4 


- x 5 


- Xl - 


h X4X5 - 


h X4X1 - 


h X5X1 + xix 2 


— X4X5X1 


— X4X1X2 - 


1- X4X5X1X2 


1 


— X4 


- x 5 


- Xl - 


1- X4X5 - 


h X4X1 - 


1- X5X1 - 


1- X1X2 


— X4X5X1 


— X4X1X2 - 


I- X4X5X1X2 


1 


- X 5 


- Xl 


- X 2 ~ 


V X5X1 - 


h x 5 x 2 - 


(- xix 2 - 


Hx 2 x 3 


- X5X1X2 


- x 5 x 2 x 3 - 


1- X5X1X2X3 


1 


- x 5 


- Xl 


- X 2 - 


V X5X1 - 


h X5X2 - 


1- xix 2 - 


[-X2X3 


— X5X1X2 


- X5X2X3 - 


1- X5X1X2X3 


1 




- x 2 


- x 3 ' 


V xix 2 - 


I-X1X3 - 


(- x 2 x 3 + X3X4 


- xix 2 x 3 


— X!X 3 X 4 - 


f X!X2X3X 4 


1 


- Xi 


- x 2 


- X3 - 


\- xix 2 - 


h X1X3 - 


1- x 2 x 3 - 


1- X3X4 


- X1X2X3 


— X1X3X4 - 


f X1X2X3X4 



° 1 - X2 — X 3 — Xi + X2X3 + X2X4 + X3X4 + X4X5 — X2X3X4 — X2X4X5 + X2X3X4X5 
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From our general formula given in the Appendix three invariants of the 5dLV map are 

Hi = xix 2 + X2X3 + X3X4 + x 4 x 5 + X5X1 - xi — x 2 - x 3 — x 4 - x 5 , 

Hi = X1X3 + X 2 X 4 + X3X5 + X4X1 + X 5 X 2 - XlX 2 X 3 — X2X3X4 — X3X4X5 

— X4X5X1 — X 5 XiX2 — X1X2X4 — X1X3X4 — X1X3X5 — X2X3X5 — X2X4X5 

+X2X3X4X5 + X3X4X5X1 + x 4 x 5 xix 2 + X5X1X2X3 + X1X2X3X4, (26) 
r = X1X2X3X4X5. 

If we form two particular combinations H 2 + 3Hi + 5 and r + Hi + 2 from these invariants, we see that the Grobncr 

(2) 

basis of these functions generates the 3rd elimination ideal of the functions {X- ' — xj}. Therefore the invariant 
variety of periodic points of period 2 is 

w«72» = {x\ H 2 + 3Hi + 5 = 0, r + Hi + 2 = 0}. 

This is an algebraic variety of dimension 3. 



D. 3 point Toda lattice map 

As we show in the Appendix, the N point Toda map is equivalent to the 2N dimensional Lotka-Volterra map. 
Since it is much easier to manage the Toda map than the corresponding Lotka-Volterra map from the computational 
point of view, we consider here the 3 point Toda map, 

( /w , i« , , v« , ) - , it 1} , 4 t+1) M t+1 \ vt 1] M t+1) ) , 

which is defined by 

r (t+i) _ At) V 3 V\ + h h + h v i v (t+i) _ v (i 

2 i^+v^v^+i^^y 1 

.(t+i) _ (t) £1 V 2 + A -/ 2 + ^ V 2 v (t+i) _ v (i 

3 v^+i^i^+i^y 2 2 

, v , , At) v {t) v {t) v (i) T {t) T (t) 

j{t+l) At) h. V 3 +V 2 ^3 + h 1 3 v (t+l) v (t _ _ 

1 /fVf + ^^ + ^Vf 3 3 7«yW + y«y« + 

This map has four invariants, 

Ti = + /<*> + !<*> + if > + V 2 (t) + V 3 (t) , 

t 2 = j^/W + 4 4) 4 4) + j^/W + vfV 2 (t) + v 2 (t V 3 (t) + viV/" 

+/ 1 (t v 2 (t) +4 t v 3 (t) +4 t v 1 (t) , 
t 3 = jW/W/w, n = v/ 1 (t v 2 (t v 3 (t) . 





+ y 2 (t v 3 (t) + 4 t) 4* ) 






i +4 t v 1 (t) ' 






, +4 t v 1 (t) 


At) v (t) 

1 l v 2 


"•" ^1 J 2 


+ y 1 (t V 2 (t) ' 


At) v (t) 

1 l V 2 


, r(')rW 

"+" 7 1 J 2 
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(2) (3) 

Prom these data we find the Grobner bases generating the 4th elimination ideal of {Xj ' — Xj} and {Xj ' — xj}. 
The invariant varieties of periodic points are thus obtained 

K<72» = {x\ T 2 = 0, T 3 -T^ = 0} 

in the case of period 2, and 

w«73» = {x\ T x = 0, T 2 = 0} 

in the case of period 3. 



V. CONCLUDING REMARKS 



Throughout this paper we have shown many examples of invariant varieties of periodic points by studying rational 
maps which are known integrable. The problem of finding a criterion which distinguishes an integrable map from 
nonintegrable one is, however, not easy. As is known, for instance, there are some subtleties involved in defining 
integrability even when all periodic points are neutral [is! ll^ . Nevertheless the existence of invariant varieties of 
periodic points seems to play an important role in characterizing integrable maps. Before closing this paper we want 
to discuss relations of our results with those well known in the classical problems of continuous dynamical systems, 
such as Poincare sections and Painleve properties. 



It is well known that sections of a periodic orbit of an integrable system draw a curve on a plane of the phase space 
as we change initial values continuously 17j. Moreover every point on an Arnold torus is neutral in all directions 
so that orbits are independent from each other. This is strongly related to the fact that an integrable system has 
sufficiently many constants of motion. When the system is perturbed the periodic points are frozen by Birkhoff's 
fixed point theorem. The situation is apparently similar to our problem of discrete maps. We must, however, take 
into account of the difference in the meaning of periodicity in these two cases. The period of the map is the number 
of steps before it returns, while in the case of the Poincare sections it is the number of sections of an orbit before it 
hits the same point on the plane. In the map there is no notion of orbit. 

In order to compare these two objects directly we may introduce a parameter 8 which specifies the interval between 
two steps of the map. If we rescale the dynamical variables X^s properly we shall obtain various integrable differential 
equations in the limit of 5 — > During this standard procedure, however, we encounter the difficulty of 

distinguishing two points in a periodic map. Consider, for example, a map of period 3 

xw -» X<$ -> X& _» x« -» • • • 

After t — 3N5 we will find the point at X^\ but see also X^> and X^ at t + 5 and t + 26, almost the same time in 
the limit of 8 — > 0. 
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Instead of introducing 5 we may interpolate the discrete map by certain smooth functions. In doing this we must 
be very careful not to destroy integrability. If we are considering maps of the Toda family we know how to manage. 
Namely we write every dynamical variable in terms of a r function of the KP hierarchy. If we simply identify the 
number of steps n with continuous time t the interpolation is established. Assuming our view is correct we should be 
able to relate our algebraic varieties obtained in this paper with solutions of the KP hierarchy given by hyperelliptic 
functions. The famous Poincare section of the 3-point Toda lattice found in |hg| . for example, must be reproduced 
as a slice of this variety. 

The Painleve test for an ordinary 2nd order differential equation tells us that it is nonintegrable if there exists a 
branch point singularity which depends on initial values. The maps we examined in this paper contain differential 
equations whose integrability is guaranteed by this test in proper continuous limit of 6 — » |S( ■ Therefore we expect 
some correspondence to the Painleve test to remain in our approach. But, again, the procedure of taking the limit 
(5^0 obscures the meaning of the period of the map. 



Finally let us comment on other works related in some way with our results. Recently Tsuda 19] discussed the QRT 
map from a geometrical viewpoint and obtained many interesting results. In particular all periodic maps of QRT are 
classified corresponding to particular combinations of the parameters. In the paper the authors derived a number 
of difference equations whose solutions are limited to periodic functions of arbitrary initial values. These equations, 
which they called 'recurrence equations', include many higher order equations in addition to those derived from the 
QRT equation. The relations of these works with our results will be interesting to be clarified. 

A simpler model which presents the same property that we discussed in this paper was invented by S.Onozawa. 
Using that model he explored more details of the transition between integrable and nonintegrable mapsj^]. 
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Appendix: r functions, Toda lattice and Lotka-Volterra maps 
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Here we explain some relations between r-functions of the KP theory, the Toda lattice and the Lotka-Volterra 
maps. All invariants of the series of the Lotka-Volterra map will be presented explicitly. 
The t function of the KP theory satisfies the Hirota-Miwa equation, 



Tj (I + 1, m + V)Tj(l, m) — Tj(l + 1, rn)Tj(l, m + 1) — tj+i (I, m)7j_i(Z + 1, m + 1) = 0, 



(27) 



with j, l,m € Z. If we define 



Ij(l,m) := 



Tj-i(l, m)Tj(l + 1, m) 
Tj(l, m)r J -_i(/ + 1, m) ' 

they satisfy the 2 discrete time Toda lattice equations: 



vn s. _ T j+ i(Z,m)T j _i(Z,m+ 1) 



Tj(Z, m)Tj(l, m + 1) 



(28) 



/,•(/, m + 1)V,(/ + l,m) = I j+1 (l,m)Vj(l,m), 
Ij(l,m+l) + Vj-i(l + l,m) = Ij(l,m)+Vj(l,m). 

The reduction to one discrete time Toda lattice is achieved by introducing new time variable t = I + m and writing 

To obtain a finite system we simply impose conditions -fj+jv = , ^j+jv = • We call it the iV point Toda map. 
If we define X.- by 



/f = (i-xg)_ 1 )d-xg)), v® = X ®X®. x 



(29) 



they define the Lotka-Volterra map 



V 



'(1-I^)=lf(l-I J 



(30) 



When the label j of xj*^ changes from 1 to d and = X^p is satisfied we call the map a d dimensional Lotka- 

Volterra map. If d is even, the Lotka-Volterra map is equivalent to the d/2 point Toda map. 



These are the results known in the literature, especially in |ll|. Integrability of these maps is obvious since their 
solutions can be written explicitly in terms of the r functions. In \v\ a simple method is given to find invariants of 
the N point Toda map. We extend this method, in the following, to the case of the d dimensional Lotka-Volterra 
map and give invariants explicitly. To this end we first define two matrices 



R(t) 



,1-X 





V i 



(*) 



i 







1 -X, 



(I) 



1 - X 



(t) 

d-1 



i-xf) 



1G 



L(t) = 



/ n 


n 

u 




n 

u 


A d 


1 \ 


l 
















1 












U 


X 3 


-1 

1 


U 






V o 







A d-1 


1 


o / 



A straightforward calculation will show that the Lotka-Volterra equations 1)30(1 are equivalent to the following matrix 
formula, 



To find invariants we define 



A{t) := L(t)R(t) = 



L(t+l)R(t + l) = R{t)L{t). 



(I 

1 1 

P2 1 1 

p 3 1 1 



Pd 1 \ 

pi 






Vo 



1 1 
p d -i 1 1 / 



where we used the notation pj :— X- (1 — X^\). Since it satisfies 

A(t + 1) = L(t + l)R(t + 1) = 
eigenvalues of A are invariant. If we write 

d 

det(A-X) = (-l) d - 1 Y / H k (X-l) k , 
the set of coefficients H$, Hi, Hd are also invariant. Comparing the both sides we find 

1- {-l) d PiP2---Pd, k = 
. . .PhPh-"P3h> k = l,2,...,[d/2] 

31,32, — ,3k 



0. 



k = [d/2] + l,...,d- 1 
k = d. 



(31) 



(32) 



Here [d/2] = d/2 if d is even and [d/2] = (d — l)/2 if (i is odd. The prime in the summation Y^, of 1(32(1 means that 
the summation must be taken over all possible combinations of ji, ji, ■■-,jk but not including direct neighbours. For 
example, in the case of d = 5, 



H 2 = P1P3 + P2P4 + P3P5 + PiPl + PbP2, 
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from which we obtain (|26|l . The number of independent invariants of d dimensional Lotka-Volterra map is [d/2] + 1. 
Since Hq can be represented by other Hk's and 

r = X i f ) X$ ) ---Xf, 

it is convenient to use r instead of Hq. 
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